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Fixed-energy harmonic functions 

Aaron Abrams * Richard Kenyon^ 
March 14, 2017 


Abstract 

We study the map from conductances to edge energies for harmonic functions on 
finite graphs with Dirichlet boundary conditions. We prove that for any compatible 
acyclic orientation and choice of energies there is a unique choice of conductances 
such that the associated harmonic function realizes those orientations and energies. 

We call the associated function enharmonic. For rational energies and boundary data 
the Galois group of (the totally real algebraic numbers) over Q permutes the 
enharmonic functions, acting on the set of compatible acyclic orientations. 

For planar graphs there is an enharmonic conjugate function; together these form 
the real and imaginary parts of a “fixed energy” analytic function. In the planar 
scaling limit for (and the fixed south/west orientation), these functions satisfy a 
nonlinear analog of the Cauchy-Riemann equations, namely 

UxVy = 1 

— 1 * 

We give an analog of the Riemann mapping theorem for these functions, as well as a 
variational approach to finding solutions in both the discrete and continuous settings. 

1 Overview 

1.1 The Dirichlet problem 

The classical Dirichlet problem is to find a harmonic function on a domain that takes 
specified values on the boundary. It is the mathematical abstraction of physical equi- 
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librium problems arising in electromagnetism, fluid flow, gravitation, and other areas. 
Solutions to the Dirichlet problem describe physical phenomena from the vibrations 
of a drum to the dissipation of heat. Dirichlet himself was investigating the stability 
of the solar system In the century and a half since Dirichlet’s work, there have 
been countless investigations into various versions, modifications, generalizations, and 
special cases of the Dirichlet problem. 

We study the discrete Dirichlet problem on a graph, and in particular we are 
interested in the connection between the edge conductanees and the edge energies 
of the resulting solution. In the classical problem one fixes the conductances, and 
then minimizing the Dirichlet energy gives a unique solution to the Dirichlet problem, 
namely a harmonic function on the vertices called the voltage. The edge energies are 
then obtained from this solution. Thus there is a map, soon called T, from edge 
conductances to edge energies. 

Here we reverse the problem and fix the desired energies of the harmonic function, 
determining which choices of conductances lead to these energies. We show (Theorem 

that in suitable coordinates T has a surprisingly simple Jacobian, which makes 
analysis of the point preimages very clean. This set of preimages is parameterized by 
a certain set of acyclic orientations of the graph (Theorems and [^, and one can 
determine the Dirichlet solutions through a simple equation, the enharmonie Laplace 
equation (Theorem]^. Remarkably, the cardinality of the solution set is independent 
of the boundary values, as long as these are distinct (Corollary]^. 

We begin in Section 1.2 with notation and an overview of our main results showing 
existence of solutions and giving variational descriptions of these solutions. In Sections 


1.3, 1.4 1.5 1.6 we give consequences of our theorems that highlight connections 


between the Dirichlet problem on a graph and other areas of mathematics, such as 
number theory (Theorem , discrete analytic functions (Theorem and rectangle 
tilings (Corollaries and [^. Some proofs are given; longer ones are deferred to later 
sections. 

Several explicit examples of the main theorems are given in Section 
In Section we study the Jacobian of 'k. Sections EEE contain proofs. We pose 
a few questions in Section 


1.2 Main results 

Let Q = {V, E) be a connected graph and B C V a subset of > 2 vertices, which we 
refer to as boundary vertices. We assume that each edge of Q lies on a simple path 
from a vertex of B to another vertex of B. Let Ce > 0 be a conductance for each edge 
e. It will sometimes be convenient to refer to edges by their endpoints, e.g., we also 
write Cxy for the same conductance if e is the edge with endpoints x, y. 

Given a function u : R —)• M the classical Dirichlet problem is to find a function / 
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equal to on i? and which minimizes the Dirichlet energy 

S{f)= Yj Ce(/(a;) -/(2/))^ 

e=xy 

where the sum is over all edges of G. There is a unique minimizing function h; it is 
the unique function harmonic at each vertex of V \B, that is satisfying, for x ^ B, 

0 = Ah{x) := Y Ce{h{x) - h{y)) 

y^x 

where the sum is over neighbors y of x. We define the energy of h on an edge e = xy 
to be Ce{h{x) — h{y))‘^, and 

: (0, oo)® [0, oo)^ (1) 

to be the map from the set of edge conductances to the set of energies of the associated 
harmonic function h. See an example in Q below. 

Our main theorem relates the conductances and edge energies of a harmonic func¬ 
tion. (See the definition of compatible orientation below.) 

Theorem 1 (Existence). Let the graph Q, boundary B, and boundary values u : i? —)• M 
be fixed. For each compatible orientation a of G and tuple £ G (0, oo)® there is a unique 
choice of positive conductances {cg} on edges so that the associated harmonic function 
h has sgn{dh) = a and energies £. 

Let Tu be the set of real-valued functions on G which have no interior extrema, 
that is, no local maxima or local minima except on B. For any choice of positive 
conductances, a harmonic function with boundary values u on B necessarily lies in 
J-u, by the maximum principle for harmonic functions. 

Suppose that f € iFu takes distinct values at the two endpoints of each edge of G- 
Then / induces an orientation of the edges of G, from larger values to smaller values. 
This orientation has the following three obvious properties: 

1. it is acyclic, that is, there are no oriented cycles, 

2. it has no interior sinks or sources (as these would correspond to interior extrema 
of /), 

3. there are no oriented paths from lower boundary values to higher (or equal) 
boundary values. 

Given u, we call an orientation satisfying the above properties compatible (with u), 
and let be the set of orientations compatible with u. It is not hard to see that for 
each fj G Su there is a function h ^ Fu inducing cr, that is, with dh{e) := h{x) — h{y) 
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nonzero on every edge e = xy, and having sign sgn((i/i(e)) = a{e) on each edge. The 
set Fu{(y) of functions h G Tu inducing the orientation a is therefore an open polytope 
since it is defined by (strict) linear inequalities. 

The most basic setting of Theorem above is the case where the boundary B 
consists in just two vertices. Let B = {uo,^!} and u{vo) = 0,n(ui) = 1. Then 
consists of acyclic orientations which have a unique sink, at uq, and a unique source, at 
ui. These are called bipolar orientations. Their cardinality is the beta invariant of the 
graph G, in which an edge is added to Q from vq to vi if there is no such edge already 
present. The beta invariant, also known as the chromatic invariant, is the derivative 
at 1 of the chromatic polynomial of or equivalently the coefficient of x (and of y) 
in the Tutte polynomial T{x,y), see [Ij. 

By Theorem above, if we fix a graph Q with boundary B, a function tt : B —)■ M, 
and on each edge e an energy £e > 0, the set of functions in Tu with energies {Sg} is 
in bijection with the set of compatible orientations of G. We call these functions 
enharmonic (short for “energy-harmonic”) with respect to the energies {£e}- Note 
that a function that is enharmonic with respect to certain energies is also harmonic 
with respect to certain conductances, and vice versa; the difference in terminology 
indicates only which data were used to define the function. 

Theorem 2 (Characterization). Let the graph G, boundary B, and boundary values 
u : B be fixed. Enharmonic functions with energies £ G (0, oo)® are precisely the 
local maxima on Tu of the functional 


M(h)^ n (2) 

e=xyGE 


The enharmonic function h with energies £ and satisfying dh = a is the unique maxi¬ 
mizer of M{h) on the polytope J-u{o'). An enharmonic function is one which satisfies 
the (discrete) enharmonic Laplace equation Lh{x) = 0 at every interior vertex x, where 
L is the nonlinear operator 


Lh{x) 


E 

y^x 


£ 


xy 


h{x) - h{y) 


(3) 


and the sum is taken over all vertices y adjacent to x. 


In practice, for given energies {£e}, determining the function h and the conduc¬ 
tances {ce} is most easily done via the equation Q. We prove Theorem]^ in Section 

El 


1.3 The degree of ^ 

The map of ([^ from conductances to energies is a rational map (i.e. its coordinate 
functions are ratios of polynomial functions of the conductances). To see this, note first 
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that to obtain the harmonic function f as a function of the conductances and boundary 
data, we have to invert the Laplacian which has entries which are linear polynomials 
in the conductances; the energies are then obtained from Se = Ce{f{x) — /(y))^. Thus 
can be extended to a rational map 

The inverse map is algebraic. Remarkably, for positive real energies all inverse 
images are positive real: 

Theorem 3. The rational map is of degree |S„| over C. 

By this we mean that there is an open dense set of C® such that the cardinality of 
the preimage of any point in this set is | |. See Section for a proof. Since Theorem 

[^provides |S„| preimages of any point with positive real coordinates, such a preimage 
must consist only of points of with positive real coordinates: 

^«^((0,oo)^) C (0,oo)^. 


An immediate consequence is 

Corollary 4. The number of compatible orientations |Su| does not depend on the 
values of u, as long as u takes distinct values at distinct boundary vertices. 

Proof. The degree of a rational family of rational maps is constant for almost every 
parameter value. Thus for almost every choice of values of u, the degree of is 
constant. As u varies among a (full-dimensional) set of functions taking distinct values, 
the set of compatible orientations does not change. So the set of u’s where the degree 
of Tu is smaller is a subset of the set u’s where at least two boundary values are the 
same. □ 

1.4 Galois action 

Suppose the given energies and boundary values lie in Q, with positive energies. Then 
the enharmonic equation 

V- — -= 0 (4) 

2. hix) - h{y) ^ ^ 

supplemented with the boundary equations h\B = u which define the enharmonic 
functions, form a rational system of equations defined over Q. Thus the absolute 
Galois group permutes the solutions. Theorem shows that the values of / generate 
a totally real number field. Thus we have 
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Corollary 5. If the values of u and the coordinates of S £ ( 0 , 00 )® are rational 
then the Galois group of Q*’'' (the totally real algebraic numbers) over Q permutes the 
enharmonic functions of Theorem 0 Eguivalently, Gal(Q*^/Q) acts on 

We conjecture that for 3-connected graphs, the action is transitiv^ 

Conjecture 1. If G is 3-connected the above Galois action is transitive for generi^ 
rational energies and boundary data. 


Considering all acyclic orientations as G and u vary, this is reminiscent of 
dessins d’enfants, a family of combinatorial objects on which the absolute Galois group 
Gal(Q/Q) acts. In the current setup it is the Galois group of the totally real algebraic 
numbers which acts. We prove (see Section that the action is faithful, i.e., that 
every totally real number field arises in this way: 


Theorem 6. For every totally real number field K there is a tuple {Q, B,u,{Se}) 
consisting of a graph Q, boundary B, rational boundary values u and rational edge 
energies Se > 0 such that K is the number field generated by the values of one of the 
associated enharmonic functions, and all enharmonic functions generate isomorphic 
number fields. 


We don’t know if Theorem holds when the size of the boundary is constrained, 
for example when the boundary consists in only two points. We show this however for 


quadratic fields in Section 2.2 


1.5 Discrete analytic functions 

A circular planar network [5] is a finite embedded planar graph Q with conductances 
c : E ^ ]R>o and boundary B consisting of a subset of vertices on the outer face. We 
define the dual network Q* as the usual planar dual of the graph obtained from Q by 
adding disjoint rays from each boundary vertex of Q to infinity. Thus Q* has one vertex 
in every bounded face of Q and one vertex in the exterior face between every pair of 
“consecutive” boundary vertices of Q ; there is a dual edge connecting dual vertices for 
every primal edge separating the associated faces and there is a dual edge for each of 
the infinite rays. Given a harmonic function f on G with boundary values u on B, a 
function g on G* is said to be a harmonic conjugate if for any primal edge e = xy we 
have f{y) — f{x) = Ce{g{t) — g{s)) where where s,t are the faces adjacent to edge xy, 
with t being on the right when the edge is traversed from x to y. These equations are 

^ A graph with boundary is 3-connected if, upon adding edges connecting every pair of boundary vertices, 
it cannot be disconnected by removing two vertices. 

^By generic we mean on a Zariski open set. 


6 




the discrete Cauchy-Riemann equations. They imply in fact that / is harmonic on Q, 
and g is harmonic on the dual graph with reciprocal conductances on the dual edges. 
In this setting such a pair /, g is also said to form a discrete analytic function. 

For a circular planar network {0,B,u,{£e}) with enharmonic function /, an en¬ 
harmonic conjugate g is a, function g on Q* satisfying 


{g{t) - g{s)){f{y) - f{x)) = S^y, 


where s, t are the faces adjacent to edge xy, with t being on the right when the 
edge is traversed from x to y. Enharmonicity of / implies enharmonicity of g on the 
dual graph with the same energies. Together / and g are harmonic conjugates for 
the associated conductances. We call the pair {f,g) the fixed-energy discrete analytic 
function associated to the boundary data. These functions are used as x- and y- 
coordinates in the rectangle tilings discussed in the next section. 


It is natural to consider scaling limits of enharmonic functions, that is, continuous 
limits of enharmonic functions on a sequence of graphs approximating (in an appropri¬ 
ate sense) a region in M"'. These scaling limits depend strongly on both the underlying 
graph structure and choice of orientation. We consider here only the simplest nontriv¬ 
ial case where the approximating graph is a scaled subgraph of , with edges oriented 
south and west. We prove in Theorem 11 below that the natural scaling limit of an 
enharmonic function / and its enharmonic conjugate g for Q = el? satisfy a nonlinear 
analog of the Cauchy-Riemann equation: 


fx9y = -1 (5) 

fy9x — 1- 


Moreover /, g both satisfy the “continuous enharmonic Laplace equation” 


fxx fyy 

f2 f2 

J X J y 


We prove (see Section 5.3) an analog of the Riemann mapping theorem in this context, 
namely that from any region in an appropriate class there is a “fixed-energy analytic” 
map to a rectangle. 

In [T] a different continuous analog of the hxed-energy problem is considered. 


1.6 Tilings and networks 

In 1940, Brooks, Smith, Stone, and Tutte [3], building on work of Dehn m , gave a 
correspondence between the Dirichlet problem for circular planar networks (with two 
boundary vertices) and tilings of a rectangle by rectangles; they called this the Smith 


7 




diagram of the network. See Figure below for an example. Namely, given such a 
tiling, one can orient the rectangle so it sits on the x-axis, and then form a network 
from the tiling according to the following rules. Place a vertex at height y for each 
maximal horizontal segment in the tiling at height y; then, for each rectangular tile add 
an edge to the network connecting the vertices corresponding to the top and bottom 
of the tile. On this edge place a conductance which is the aspect ratio (width over 
height) of the tile. The result is a circular planar graph which can be interpreted as a 
resistor network. 

If we now imagine the top and bottom vertices connected by a battery generating 
a voltage differential equal to the height of the rectangle, then the voltage function on 
the vertices (i.e., the y-values) is a harmonic function that solves the Dirichlet problem 
on the network with boundary conditions enforced at the top and bottom vertex. The 
rectangle widths are the current flows Cedh, and the areas are the edge energies Ce{dh)'^. 

The BSST correspondence is invertible: beginning with a planar network with a 
specified source and sink on the outer face and a specified total voltage drop, one can 
solve for the harmonic function and the currents in the wires. The resulting network 
can then be represented as a rectangle tiling by reversing the process described above. 

A corollary of Theorem applies to Smith diagrams. This was originally proved in 
a different manner in [12] and a number of times since, see |8|. 

Corollary 7. Given a rectangle tiling of a rectangle there is an isotopic tiling in which 
the rectangles have prescribed areas. 

Here by isotopic we mean one is obtained from the other by displacing continuously 
the walls, floors and ceilings of rectangles in the first to form the second, without 
degenerating any rectangles along the way (we also need to add a proviso that at a 
point where four rectangles meet only one of the two possible perturbations is allowed; 
either the horizontal segments to its left and right can move independently up and 
down, or the vertical segments above and below can move independently left and 
right, but not both.) See Section]^ for some examples. 

Such maps are also called rectangular cartograms |9] and floorplans |2] and are used 
e.g. to make geographic maps where the area represents the measure of some quantity 
associated to each country, such as population. 

Another corollary of Theorems and applies to the tilability of polygons by 
rectangles of rational areas. Such a tiling arises from a circular planar network with 
a boundary vertex for each horizontal segment of the boundary polygon. Again the 
Dirichlet boundary values are the y-coordinates of the horizontal boundary segments. 

Corollary 8. If a polygon P can be tiled with rectangles of rational areas, then the 
horizontal edge lengths are in a totally real extension field of Q{vi,... ,Vk), the field 
generated by the vertical edge lengths. 
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Figure 1: A polygon with area 1 which is not tilable with rational-area rectangles. 

For example the unit-area polygon of Figure cannot be tiled with rectangles of 
rational areas, because 2^/^ is not a totally real algebraic number. 

Proof. Assume some vertex of P is at (0,0). As discussed above, the y-coordinates 
of the horizontal edges of P correspond to the boundary values u of the enharmonic 
function. The y-coordinates of the interior horizontal edges are the interior values of 
the enharmonic function, and therefore in a totally real extension of Q(ui,..., u^). 
The x-coordinates of all vertical edges are sums of currents £xy/{f{x) — f{y)), which 
are therefore in the same extension field. □ 


2 Examples 

We illustrate our theorems and techniques with a few examples. 

2.1 A small graph 

Consider the graph G shown in Figure [2| with conductances a, b, c, d, e as labeled. 
With u{vi) = 1 and u{vo) = 0, we have x,y G [0,1], and thus Tu — [0,1]^. Solving 
for h gives h{x) = ^{ab + ac + be + ae) and h{y) = ^{ab + ac + be + bd), where 
Z = ab + ae + ae + be + bd + ed + ee + de IS the determinant of the Laplacian. The 
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Figure 2: A graph with |S| = 2. With u{vo) = 0,m(ui) = 1 and all energies equal, the 
underlying number held is Q(\/5). 


energies are then 

b, c, d, e) = -;^(a{bd + cd + ce + de)^, b{ac + cd + ce + de)^, 

^ '■ ( 6 ) 

c{ae — bd)‘^, d{ab + ac + be + ae)'^, e{ab + ac + be + bd)"^ ] 


To find solutions with (for instance) all energies equal to 1, it is much easier to 
solve the enharmonic equation Lh = 0 than it is to set each coordinate above equal to 
1 and solve the system directly. The enharmonic equation(s) (see Q) are 


£a £c £d _ 1 1 1 

h{x) — h{vi) h{x) — h{y) h{x) — h{vo) h{x) — 1 h{x) — h{y) h{x) 

£b Sc Se 1 1 1 

h{y) - h{vi) h{y) - h{x) h{y) - h{vo) h{y) - 1 h{y) - h{x) h{y) 


and have the solutions h{x) = ^ ± ^ and h{y) = 5 T from which one can quickly 
produce the conductances (for example, a = e = ^(15 ± 5 \/ 5 )). The values of h are 
roots of the polynomial 5z‘^ — 5^; + 1. There are two solutions because current can flow 
either way along the edge xy. That is, |S| = 2. The two polytopes in J^u = [0,1]^ on 
which M(h) is nonzero are the triangles above and below the diagonal. 

The equal-area tilings associated to the two solutions are shown in Figure the 
y-coordinates of the horizontal edges are h{x) and h{y). The tilings are not isotopic. 


2.2 Quadratic number fields 

Keeping the same graph Q as in the previous example, if we instead look for solutions 
with arbitrary energies £a,Sb,Sc,Sd,Se, we find that the conductances generate the 


10 
























number field Q(\/5), where 


5 = {8a£e - £b£d + £cSf + ^£b£c£dS. 

Here S is the sum of the five energies. Of course for positive real energies, 5 > 0 as 
all roots are real. It turns out that we can get any real quadratic number field this 
way. To see this, we first specify £a = £b = £c = ^ so that the discriminant reduces to 
5 = dTj+dT^Te+dTe+12T(i+12<£’e+9. To obtain the number field for a positive 

squarefree integer D, it suffices to find positive rational numbers £d, £e such that 
S = Dr^ for some rational r. Choose a rational s so that D' = Ds^ G (1/3,4/9). Then 
substituting £d = = \~iy -2 (both positive) gives discriminant 5 = 

which is a rational square times D. 

2.3 Jacobi polynomials 

One way to generalize the previous graph is the followingj^ Let Q = Kn +2 be the 
complete graph on n + 2 vertices uq, ... ,u„+i, minus the edge vqVi. Fix boundary 
B = {uo,ui} and boundary values m(uo) = 0 and u{vi) = 1. Put energy a on all edges 
from vq to Vi, z G [2, n + 1] and energy b on edges from vi to Vi, z G [2,n + 1], and 
energy 2 on the remaining edges. Then by symmetry we can order the values h{i) for 
z G [2, n + 1] in increasing order. We have 

M{h)= n Hifd - Hi))" n (Hi)-Hi))"- 

i&[2,n+l] l<i<j<n+l 

In m Theorem 6.7.1] it is shown that the values h{i) maximizing M{h) are the roots 
of the Jacobi orthogonal polynomial scaled to the interval [0, Ij. 


2.4 A medium graph 

For another example, consider the network shown in Figure]^ Again let u{vi) = 1 and 
u{vo) = 0. The edges are labeled 1 through 12, and (by coincidence) also |S„| = 12. We 
have solved the enharmonic equation (see (Q) for the twelve enharmonic functions that 
produce all energies equal to 1. From these functions we computed the conductances of 
the networks, which determine the shapes of the rectangles in the corresponding Smith 
diagrams. In Figure]^ we have scaled the Smith diagrams to be (unit) squares, so each 
rectangle in each rectangulation has area 1/12. Any rectangulation of a rectangle with 
the same underlying network (Figure]^ is isotopic to one of these twelve. 

^Note that n = 2 is the only member of this family that is planar. 
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These calculations are too complicated to do by hand; for instance, according to 
Mathematica, the width of the rectangle labeled “1” in each picture is a root of the 
polynomial 

p{z) = 1270080000000^^2 - 5554584000000^^^ + 10776143400000z^° 

- 12235337185000^^ + 9034493949125z® - 4560532680000z^ 

+ 1610724560815Z® - 400501165895^® + 69535433439^^ 

- 8223166134Z® + 629396649^2 - 28041714z + 551124. 

All other edge lengths in all tilings in Figure lie in the splitting field of p, which is 
contained in M. 



V 


Figure 3: A planar network with two boundary vertices no,ni. 


2.5 A large graph 

The variational method makes it very easy to find (numerically) enharmonic functions 
on large graphs, see for example FigureIn fact this figure illustrates a certain scaling 
limit convergence theorem. Theorem |11| below. 

3 The Jacobian 

The work in this paper began with the observation that the Jacobian of T has a 
surprisingly simple description, which we give here. Our proof of Theorem uses the 
nonvanishing of this Jacobian. 
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Figure 4: Smith diagrams for the planar network of Figure]^ with energies 1. Here |S„| = 12, 
and the Galois group acts transitively: the underlying number held for each tiling is a degree- 
12 totally real extension held of Q. Each tiling corresponds to a diherent bipolar orientation 
of Q, so these tilings are not isotopic. 
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Figure 5: Smith diagram for the 40 x 40 square grid [0,40]^ fl with edges oriented south 
and west, with all energies 1 and boundary values m(0,0) = 0, m( 40,40) = 1. The apparent 


regularity here can be quantihed: see Theorem 11 
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First some definitions. On a connected graph Q = {V,E) with boundary vertices 
B, a 1-form is a function to on oriented edges which is antisymmetric under changing 
orientation: uj(—e) = —co(e} where —e represents the edge e with reversed orientation. 
Let be the space of 1-forms. Let F^ be the space of functions on V with zero 
boundary values. 

We define the map d : F^ —)■ F^ by df{e) = f{y) — f{x) when e = xy is directed 
from X to y. A 1-form is a coboundary if it is of the form df for some / G F^. Let 
Wcob C F^ be the space of coboundaries. It has a natural basis consisting of the set of 
coboundaries {dlv}v&v\B where is the function which is 1 at u and zero elsewhere. 

The orthogonal complement (for the standard inner product on M^) of Wcob is 
the space Wcyc of (relative) cycles, or flows. These are 1-forms in the kernel of the 
boundary operator d* ; F^ —)• = F*^, the transpose of d. The mapping d* from 

1 -forms to functions is 

d*u}{x) = ^iw{xy) 

y^x 

where the sum is over vertices y neighboring x; in other terminology d*oj is the diver¬ 
gence of Lo. Cycles (or flows) are 1-forms which are divergence-free on the non-boundary 
vertices. 

Note that Wcob has dimension \ V\ — \B\, the number of internal vertices, and Wcyc 
has dimension \E\ — \V\ \B\. 

The differential of the map Tm has a surprisingly simple form when considered as 
a map from the logarithms of the conductances to the logarithms of the energies. We 
assume the boundary values u are fixed. If we fix an element cr G we can identify 
with the space of 1 -forms. 

Let Pcyc : F^ —)■ F^ be the orthogonal projection onto Wcyc, and Pcob = I — Pcyc the 
projection onto Wcob- 

Theorem 9. When T^dc}) is nonzero on every edge, let Jiog he the differential of 
the map logoT^ o exp : —>■ M®. Then flog = Pcyc — Pcob- 

Corollary 10. 

detUT = (-l)l^l-l^l Yl{h{x) - h{y)f. 

e 

In particular on (0, oo)®, T fails to be loeally injective preeisely when some edge has 
energy 0. 

Proof of Corollary \1C\ By the Theorem (and the chain rule) 

detH^ = (_i)d*m(W^co 6 )n^. 

rieCe 

Substitute £e = Ce{h{x) — h{y))‘^ for edges e = xy. □ 
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Proof of Theorem^^ Since = Wcob®Wcyc, we study the action of Jiog on Wcob and 
W^cyc- 

Given a conductance function c, let h be the associated harmonic function, and 
oj = c ■ dh the associated current flow. We assume cu is nonzero on each edge, which 
holds for almost all choices of c. 

Let V G V \ B. We define a 1-parameter family of conductances Ce{t), t G [0,1] as 
follows. Let hi{v) be any real number close enough to h{v) so that, if we redefine h{v) 
to be hi{v) then the sign of dh does not change. Fix h at all vertices but v and let 

ht{v) = h{v) + t{hi{v) - h{v)). 

We fix Lo to its initial value, and define the conductance at time t by Ce{t) = u{e)/dht{e). 
Then ht is harmonic for conductances Ce{t). The associated energies are £e{t) = 
uj{e)dht{e) and thus ^logTe(^) = “^logCe(t). Evaluating at t = 0, this vector is a 
multiple of dl^. 


Ft 


log Ce{t) 


t=0 


m 


log dht{e) 


t=o 


hi{v) — h{v) 
h{v) 


dly. 


Thus dly is an eigenvector of Jiog with eigenvalue —1. As u varies over V \ B the 
coboundaries dly span the coboundary space Wcob, so Jiog is minus the identity on 
Wcob. 

Now let 7 be a cycle in G. Let 5 > 0 be small enough so that ui = uj + 6'^ has the 
same sign as lo. We now fix h but vary to according to LOt{e) = uj{e) + t{uji{e) — Lo{e)). 
Define the conductance at time t by Ce(t) = uJt{e)/dh{e). We then have ^logTe = 
^ logCe(t). Evaluating at t = 0, this vector is a multiple of 7 : 


Ft 


logCe(t) 


t=0 


Ft 


loga;i(e) 


t=o 


a;i(e) — a;(e) 
uj{e) 


7- 


Thus Jiog is the identity on Wcyc- 


□ 


It is worth pointing out the interpretation of this proof in terms of tilings. If G is 
planar, the proof starts with a rectangle tiling with combinatorial structure described 
by G, and then studies the effect of isotopy on the diagram. Moving a horizontal 
segment corresponds to changing the value of /i at a single point, which is an eigen¬ 
value — I of Jiog: the relative change of the conductances (widths over heights) move 
in the opposite direction as the relative change of areas (widths times heights). Mov¬ 
ing a vertical segment corresponds to the eigenvalue -|-I: the relative change of the 
conductances move in the same direction as the relative change of areas. 
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4 Proofs of Theorems [T] and |3 


We assume u is fixed and takes distinct values at distinct boundary vertices, and we 
omit it from the notation. 

Proof of Theorem\^ We first prove that the map is surjective onto S, that is, every 
u G S is realized by some choice of conductances. Given u G S, let / be any function 
in Such a function exists because a is acyclic: for example one can define / 

from a linear extension of a, thinking of cr as a partial order. There is a choice of 
conductances for which this / is harmonic: to see this take any flow uj with orientation 
a and divergence 0 on V \ B, and define the conductance Cg on edge e = xy hy 
Ce = f{x)^-\{y) ■ ^ concrete example of such a flow w, let P be the set of all 

possible oriented paths from a boundary vertex to a boundary vertex, with orientation 
compatible with cr, and let the value uj{e) be the probability that e is on a uniformly 
chosen path in P. By our condition on Q that each edge lies on a simple path from B 
to B, we see that u is nonzero. This completes the proof that the map is surjective 
onto S. 

Now consider the map T: (0, oo)® —>■ [0, oo)^ sending a tuple of conductances to 
the corresponding tuple of energies of the associated harmonic function. Corollary [T0| 
shows that the Jacobian of 'k is nonzero when the energies are nonzero. 

We claim that T is proper: the preimage of a compact set is compact. In other 
words, 'k maps the boundary to the boundary. Take a sequence of conductance func¬ 
tions leaving every compact set of (0,oo)®. By taking a subsequence, we can assume 
it converges to a boundary point of [0, oo]^. We need to show that along this sequence 
some edge has energy tending to 0 or oo. 

Let So be the set of edges whose conductances are going to zero, and Soo those whose 
conductances are going to oo. If Soo does not connect any two boundary vertices, then 
there is a limiting network obtained by contracting all edges in Soo- The resulting 
network has conductances bounded above, which implies a finite current flow. The 
energies on the contracted edges and any edges of zero conductance have gone to zero. 

If Soo contains a connection between two boundary points (with different values of 
u, by genericity), the current flowing between these points will diverge. This implies 
that the energy (which can be written as current flow times dh) will diverge along 
some edges of this connection. 

Thus T is proper. 

The preimage of (0, oo)^ in the space of conductances is a union of open sets C/o- 
where dh has constant sign u. At a finite boundary point of any such Ua, some energy 
goes to zero (since zero current implies zero energy). Hence the boundary of maps 
into the boundary of the space of energies. Thus T on any of these open sets Ua is 
proper and a local homeomorphism (a covering map), and therefore surjective from 
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this set to the space of energies. To get uniqueness it remains to show that is 
connected. 

In fact we show that has the structure of (the interior of) a polytope. As in the 
first paragraph of the proof, a point in Ua is determined by a flow compatible with 
a (and of divergence zero on interior vertices), and an independent choice of function 
/ in the open polytope The set of current flows ui compatible with a is the 

subset of (0, oo)^ defined by the linear equations div(a;) = 0. It is the interior of a 
convex polytope (possibly inhnite but with finitely many facets). □ 

Proof of Theorem\^ We have constructed |S| real preimages of each point in (0, oo)®. 
It remains to show that there are no non-real preimages. Let {cg} £ T“^({Te}) be a 
preimage of a set of positive real energies. Let h be the corresponding solution to the 
Dirichlet problem with conductances {cg}. Then for all non-boundary vertices x, 

^Ce(/i(x) - h{y)) = 0 

y'^x 

and so 

_ n 

^ h{x) - h{y) 

We can rewrite this as _ _ 

^ Txy(h(x) - h{y)) ^ 

\Kx)-h{y)\^ 

or 

(LJ \h{x)Ph{:ti)\^ 

This represents h[x) as a convex combination of its neighboring values h{y), and thus 
h{x) lies in the convex hull of the neighboring values h{y). Since the boundary values 
of h are real, each h{x) must be real as well. Finally Cg = Tg/ (/i(x) — h{y))‘^ so each Cg 
is real and positive. □ 

5 Enharmonic functions: a variational method 

5.1 Proof of Theorem [2] 

Proof of Theorem\^ Fix a compatible orientation a. Recall the open polytope J-’uic) C 
Tu of functions inducing orientation a. Restricted to T'„((t), the function log M(/i) is 
concave, since it is a sum of concave functions 

£{e)\og\h{x) - h{y)\ = £{e)\og{a{e){h{x) - h{y))). 
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The function \ogM{h) tends to —oo on the boundary of T'u(o') (when h{x) = h{y) 
for some edge e = xy). As a consequence log M{h), being analytic, is strictly concave 
on Tu(it). It thus has a unique maximizer on the interior of J-'u(cj), for each cr. The 
maximizer is determined by setting the derivative with respect to h{x) equal to zero 
for each vertex x, that is 


y_ ^ _= 0. 

Kx) - h{y) 


(7) 


Given a solution to set Ce = Sel{h{x) — h{y))‘^, then 0 becomes 


'^Ce{h{x) - h{y)) = 0, ( 8 ) 


that is, h is harmonic for the conductances Cg and has energies Eg. 

□ 


5.2 Scaling limits 

Let T) C be a Jordan domain with piecewise smooth boundary and let u : dD —)• M 
be continuous. We assume u is feasible in the sense that the space ^l{D,u) of smooth 
functions /: D —)■ M with boundary values /{qu = u and satisfying fx,fy > 0 is 
nonemptjQ Fix such a function /o G Q{D,u). 

For sufficiently small e > 0 let = e7? (ID he the part of the graph el? contained 
in D. Let a be the south and west orientation of the edges of Ue. Let dD^ be the set 
of boundary vertices of De, that is, vertices with a nearest neighbor in el? not in 
Let Ue : dD), —)> M be the restriction of /o to dD,,. 

Let /e be the unique enharmonic function on with energy per edge, with 
boundary values Ue and with dfe having orientation a. 

Theorem 11. Under the above hypotheses, as e —>■ 0 the function f^ converges to 
a continuous function on D, real analytic in the interior of D, and satisfying the 
enharmonic Laplace equation 

fxx fyy 

f2 ' f2 
Jx Jy 

There is a unique (up to an additive constant) conjugate function g, also real analytic 
and continuous up to the boundary, satisfying 0. 

^The issue of feasibility is not completely straightforward. Such a function / necessarily has level contours 
having tangents pointing in the NW and SE quadrants. Conversely, given a foliation T of U by SE-directed 
leaves, in which the value of u on the endpoints of each leaf are the same, there is a function / whose level 
lines are the leaves of T. 
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We call / satisfying ([^ an enharmonic function. A function g satisying ([^ is an 
enharmonic conjugate. 


Proof. Assuming f £ C'^ satisfies (9), the 1-form to = ^ is closed. Let g{p) = 

Lo for some fixed po £ D. Then / and g satisfy the “fixed-energy” Cauchy Riemann 
equations ([^. The enharmonicity of g follows from ([^ as well. 

Let F £ be the (unique) function maximizing the concave functional 


M{F)= [ \og{F^Fy)dxdy. (10) 

JD 

Existence and uniqueness of F both follow from strict concavity of M (existence 
follows from the fact that the space of monotone functions with boundary values 
u is compact, and M. is upper semicontinuous; uniqueness follows from the fact that 
M{^{Fi+F 2 )) > M.{Fi)+M{F 2 ) unless Fi = F 2 . Analyticity follows from analyticity 
of the function log, see |TU] [Theorem 1.10.4 (vi)]. 

We will show that —?• F. 

The Euler-Lagrange equation for a critical point of M{F) is the enharmonic equa¬ 
tion 

F F 

I yy — n 
/?2 7?2 

X y 

Let Mg be the functional Q for the graph Dg with energies e^. Let /g be its 
maximizer. After triangulating the faces of Hg by adding diagonals, we can extend /g 
to a piecewise linear function /g on D (except near the boundary). Then logMg(/g) is 
a Riemann sum for A4(/g), so 

logMg(/g) = M(/g) + o(l). 

If we restrict F to we likewise hnd by analyticity of F that 

logMg(F) =M{F)+o{l). 


Consequently 

M(F) = logMg(F) + o(l) <logMg(/g) + o(l) =M(/g)+o(l) <M(E) + o(l). 

Therefore all these quantities are within o(l) of each other. By upper semicontinuity 
of A4, we hnd that /g converges to F. 

We prove the second statement of the theorem in the coming section. 
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5.3 Riemann mapping to a rectangle 

Let D C be a piecewise smooth domain of the following type. There are four points 
a, b, c, d on the boundary of D in counter-clockwise order, so that the boundary of D 
between a and b is monotone in the NE direction (that is, is monotone in the -|-x, +y 
directions), and similarly between b and c, c and d, d and a respectively the boundary 
is monotone in the NW, SW, and SE directions. Let u be the function which is 1 on 
dD between b and c, and zero between d and a. Let be as before the subgraph 
of e7? with vertices in D. Let a^,..., be boundary vertices so that the boundary 
of is monotone NE, NW, SW, SE between and b^, b^ and Ce, q and d^, and d^ 
and Ug. Let D[ be the dual graph of D^, with a vertex for every bounded face of 
and, instead of one exterior vertex, four exterior vertices Vab, Vbc, Vcd, Vda outside of the 
corresponding boundary edges of D^. Let a be the south/west orientation of edges of 
D,. 

Let / G Q,{D,u) be the unique function on D with boundary values 1 on be and 
0 on da, free on the other boundaries, and maximizing Ad(/); this differs from the 
previous case because u is only hxed on part of the boundary, but again existence and 
uniqueness follow from concavity and analyticity from arguments of unj. Then / will 
satisfy (§ on D, with fixed boundary values 0 on da and 1 on be respectively and along 
the free boundary ab U cd, f will satisfy the condition that ^ ~ ^ (this is in fact 

the Euler-Lagrange condition along the free boundary). 

As before the function on maximizing M will lie within o(l) of /. Using 
edge energies e^, let be the enharmonic conjugate of /e, defined on the dual graph 
Z?' and with value zero on the dual vertex Vab- Let = geived) be the value at the 
opposite dual vertex. Since the sum of the energies is the area of the image under the 
rectangle tiling, we have 

Re = e^(# edges of De) + o(l) = 2Area(D) -|- o(l). 

The above convergence argument applied to ge shows that ge converges to g, the 
maximizer of Af with boundary values A = 2Area(D) and 0 on cd and ab respectively. 

It remains to show that / and g are enharmonic conjugates. The enharmonic 
conjugate f* of / is the enharmonic function (well defined up to an additive constant) 
which is constant along ab and along cd, with free boundary values along be and da. 
These are precisely the boundary conditions satisfied by the function g, except that 
the value of the constant along cd might be different, so we have that f* = Bg for 
some constant B. The map {x,y) i—)• (/*,/) is area preserving except for a factor 2 
(the Jacobian is identically 2, by ([^) so the image of D has area A, which means the 
boundary value of f* along be must also be A. Thus f* = g. 

This proves that {fe,ge) converges to {f,g). In the limit we have a real-analytic 
mapping {g, f) from D to the rectangle [0, A] x [0,1]. □ 
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6 Realizing number fields 

Proof of Theorem^ Let ^ be a d + 1-star graph with vertices {?;, vq, ..., Vd}, with v 
being the central vertex. Take boundary B consisting of all vertices except v. Let e* be 
the energy on edge vvi and take rational boundary values Oj at Uj, with oq < • • • < a^. 
The enharmonic equation for the value x at u is 


d 


E 



= 0 . 


( 11 ) 


This equation has exactly d roots, and these are real and interlaced with the Oj, that 
is, there is one in each interval (oj, Uj+i) (this follows from Theorem above or simply 
from the fact that the roots are roots of the derivative of — aiY^). 

Conversely let p{x) be a polynomial with rational coefficients of degree d with all 
real roots, interlaced with the a* as above. Then we claim that there is a unique choice 
of rational Cj so that the enharmonic equation 0 has exactly the same roots. 

The map from the tuple of energies {ei} to the tuple of coefficients {bi} of the 
polynomial p{x) = bo + bix + ■ ■ ■ + bdx'^ is a rational linear map from to itself. 
Specifically its matrix consists of elementary functions in the a* (after removing aj_i 
in those for the ith column) illustrated for the case d = 3: 


/bo\ 


/ —010203 

—O0O2O3 

—O0O1O3 

—O0O1O2 \ 


/eo\ 

bi 


O1O2 -|- O3O2 -|- O1O3 

O0O2 -l- 0302 -l- 0903 

OoOi -|- 0301 -|- 0903 

OoOi -|- 0201 -|- 0002 


ei 

b 2 


—Ol — 02 — 03 

—Oo — 02 — 03 

—Oo — Ol — 03 

—Oo — Ol — 02 


62 

\hj 


V 1 

1 

1 

1 y 


\e 3 / 


The determinant of this matrix is a polynomial of degree ( 2 ) and has factors o* — aj 
for each i < j, so is the Vandermonde ± ~ therefore nonzero. Thus 

the map from energies to coefficients bi is invertible over Q. This completes the proof 
of the claim. 

Now let iL be a totally real number field and x G iL a primitive element (i.e. one 
for which K = Q[x]). Let p{X) be its minimal polynomial. Let ri,..., be the roots 
of p{X), which are real. Replacing x with x + m for a large enough integer m we can 
assume that all Xj > 0. Let oi,... , 0^+1 be positive rationals interlacing them. Now 
do the above construction with parameters a*. □ 

Note that as a simple corollary to this proof, by letting the e* be arbitrary posi¬ 
tive rationals, modulo a global scale, we can parameterize the polynomials with roots 
interlaced with the a, as a simplex 


{(eo, ■■■,ed) \ = 1}. 
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7 Other questions 

1. Can one use the volume-preserving character of the map log o'!/ oexp to compute 
the cardinality of S? This is NP hard in general but integrating the Jacobian 
over projective space one can approximate the degree of T and therefore |S|. 

2. The minimal polynomials for geometric quantities (such as values of the harmonic 
function) seem to have coefficients which are (fixed-sign) integer-coefficient poly¬ 
nomials in the edge energies. Can one give a combinatorial description of these 
polynomials? 

3. What can be said about solutions to the enharmonic equation in the case when 
the boundary values u are complex? The number of solutions is the same, but 
there is no longer the same interpretation in terms of orientations. 
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